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The macroscopic quantum theory of the electromagnetic field in a dielectric medium interacting 
with a dense collection of embedded two-level atoms fails to reproduce a result that is obtained 
from an application of the classical Lorentz local-field condition. Specifically, macroscopic quantum 
electrodynamics predicts that the Lorentz redshift of the resonance frequency of the atoms will be 
enhanced by a factor of the refractive index n of the host medium. Ifowever, an enhancement factor 
of (n^ -I- 2)/3 is derived using the Bloembergen procedure in which the classical Lorentz local-field 
condition is applied to the optical Bloch equations. Both derivations are short and uncomplicated 
and are based on well-established physical theories, yet lead to contradictory results. Microscopic 
quantum electrodynamics confirms the classical local-field-based results. Then the application of 
macroscopic quantum electrodynamic theory to embedded atoms is proved false by a specific example 
in which both the correspondence principle and microscopic theory of quantum electrodynamics are 
violated. 

PACS numbers: 42.50.Nm,42.50.Ct,03.50.De 
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I. INTRODUCTION 

The effect of a dielectric host medium on the spon- 
taneous emission rate of a two-level atom remains an 
interesting and challenging problem in quantum optics 
with an importance that befits a rigorous test of our un- 
derstanding of the interaction of light with matter. The 
essential characteristic of the problem is its multi-scale 
nature with the atom being a creature of microscopic 
quantum electrodynamics and the dielectric manifesting 
in the realm of classical continuum electrodynamics. In 
principle, it is possible to represent both the atom and 
dielectric microscopically although, in practice, the spon- 
taneous emission rate is calculated using a macroscopic 
quantum electrodynamic theory [ll, [Sj, [j, H, 0| in which 
the continuous dielectric medium is incorporated into 
a medium-assisted electromagnetic field. While the di- 
electric renormalization of the spontaneous emission rate 
of an embedded atom generates considerable interest, a 
mere handful of papers discuss the effect of a dielectric on 
the Lorentz redshift of the resonance frequency of two- 
level atoms. Knoester and Mukamel [7*1 used a Hopfield 
Q polariton model of macroscopic quantum electrody- 
namics and found that the Lorentz redshift is enhanced 
by a factor of the refractive index n compared to the 
vacuum value. Crenshaw and Bowden 9] derived the en- 
hancement factor {n^ + 2)/3 using the Bloembergen 
procedure in which the classical Lorentz local-field con- 
dition is applied to the optical Bloch equations. The two 
procedures also produce different values for the renor- 
malization of the field (Rabi frequency) that drives the 
dynamics of atoms in a dielectric. 

In this paper, we derive the generalized optical Bloch 
equations for a dense collection of two-level atoms in a 
dielectric host material using, first, the Bloembergen pro- 
cedure based on the classical Lorentz local-field condition 
and, second, Ginzburg macroscopic quantum electrody- 



namics. The derivations are simple and direct and based 
on well-established physical theories, yet produce con- 
tradictory results for (i) the Lorentz redshift of the reso- 
nance frequency of two-level atoms densely embedded in 
the dielectric and (ii) the renormalization of the field that 
drives the dynamics of the embedded atoms. The micro- 
scopic quantum electrodynamic procedure described by 
Crenshaw and Bowden [ll| confirms the Bloembergen- 
based results. Then, we must conclude that the macro- 
scopic quantum electrodynamic theory applied to dielec- 
trically embedded atoms violates both the correspon- 
dence principle and microscopic quantum electrodynam- 



II. DIELECTRIC LOCAL-FIELD EFFECTS 

In the view of Lorentz, classical continuum electrody- 
namics is better expressed in terms of an atomistic model 
of discrete particles embedded in the vacuum that in- 
teract with the microscopic electromagnetic field at the 
point of the particle [H, . 



The local field that acts on 



that particle, 
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(2.1) 



is comprised of the macroscopic Maxwell field E and the 
reaction field of all other particles, expressed in terms of 
the macroscopic polarization P. For a linearly polariz- 
able material, the polarization P = pN^FiL is the prod- 
uct of the microscopic polarizability p, the dipole number 
density Nd, and the local field. Using the Lorentz local- 
field condition (|2.ip to eliminate the microscopic local 
field produces the Clausius-Mossotti-Lorentz-Lorenz re- 
lation 
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between the polarizability and the macroscopic dielectric 
constant e = 1 + AttP/E. 

The local-field principle also applies to nonlinear me- 
dia. Bloembergen [lo| investigated nonlinear optics in 
the presence of a linear host medium and found that a 
local-field dielectric enhancement factor of 

(2.3) 



accompanies each appearance of a macroscopic field in 
the nonlinear susceptibility. In our notation, e = 
for a dielectric and H. refers to the specific quantity in 
the preceding equation and not to any other constant 
or variable representing a local-field factor. Bowden and 
co-workers [ij, IH, 13 predicted intrinsic optical bista- 
bility in a dense collection of vacuum-embedded two-level 
atoms due to an inversion-dependent local-field shift of 
the resonance frequency. Later work [9[ reported the ef- 
fect of embedding the dense two-level systems in a host 
dielectric. The dynamics of the two-level systems are 
described by the optical Bloch equations [tT]. Building 
on Bloembergen's work JjOj, Bowden and Dowling [185 
showed that the field that drives the atoms is the local 
field (|2.ip . With that substitution, one obtains the gen- 
eralized optical Bloch equations [18], 
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Here, fields are represented in the plane- wave limit by 
envelope functions such that P = ^{Pe^'^'^ 



E = i(£:e- 



c.c), and E^ = \{£Le 



-H c.c.j, 
:.). The 

macroscopic spatially averaged atomic variables in a ro- 
tating frame of reference are i?2i — (P2ie*"*)sp, R12 = 
(/9i2e~*'^')sp, and W ^ R22 - Rii = {p22)s-p - (pii)sp, 
where ( • ■ • )sp corresponds to a spatial average over a vol- 
ume of the order of a resonance wavelength cubed and the 
Pij are the density matrix elements for a two-level sys- 
tem with a lower state |1} and an upper state |2}. Also, 
/i is the matrix element of the transition dipole moment, 
assumed real, is a phenomenological dipole dephas- 
ing rate, 7|| is a phenomenological population relaxation 
rate, and Woq is the population difference at equilibrium. 

For a linearly polarizable material, the polarization is 
P = pNdEiL. For atoms embedded in a linearly polar- 
izable material, the polarization is the sum of the linear 
and nonlinear components. Substituting the local field 
(12. ip into the linear component, we have the polarization 
envelope 



P ^pNd\£ + ^P 



where A'' is the number density of atoms. Collecting 
terms in P and using the Clausius-Mossotti-Lorentz- 
Lorenz relation (|2.2p yields 



P 
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Substituting the polarization envelope (j2.5p into the gen- 
eralized Bloch equations (|2.4p produces 
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For now, we assume that £ is real. Then, the local-field ef- 
fect of the dielectric is simply an enhancement of the driv- 
ing field £, or Rabi frequency fi£/h, and the inversion- 
dependent Lorcntz redshift by £ 0]. The decay rates 
remain phenomenological. 

The optical Bloch equations [13] are the basic semiclas- 
sical equations of motion for an isolated two-level system 
in the vacuum. The two-level system interacts with its 
environment through the local field. Using the classical 
Lorentz local-field condition shows the effects of the envi- 
ronment to be: (i) a Bowden-Lorentz redshift of the reso- 
nance frequency by {Att / 3fi) N fi'^W due to nearby atoms, 
(ii) a Bloembergen enhancement of the field £ by i due 
to the linear host, and (iii) a Bloembergen- type enhance- 
ment of the Bowden-Lorentz redshift by £ due to the lin- 
ear host. Further consideration of the generalized Bloch 
equations, derived using the classical Lorentz local-field 
condition, is suspended until quantum electrodynamic 
equations of motion have been derived. At that time 
Eq. (|2.6p will serve as the basis for quantum-classical 
correspondence. 



III. MACROSCOPIC QUANTUM 
ELECTRODYNAMICS 



Quantum electrodynamics can be viewed as the quan- 
tized version of Lorentzian electrodynamics in which dis- 
crete quantum particles interact with the vacuum field 
modes. When apphed to dielectrics, the practice has 
been to create a macroscopic version of quantum electro- 
dynamics along the lines of a quantized version of con- 
tinuum electrodynamics. The macroscopic theory can 
be derived either by quantizing the classical Maxwell 
fields or by applying a continuum approximation to 
the microscopic quantum electrodynamic Hamiltonian. 
Ginzburg ^ pioneered the procedure of canonical quan- 
tization of the field in a dielectric and applied it to 
Cherenkov radiation. The macroscopic quantization pro- 
cedure was limited to dielectrics with negligible disper- 
sion and absorption. Jauch and Watson |2l.ll9l] continued 
Ginzburg's work and extended macroscopic quantization 
to dispersive dielectrics [1, [13], while other researchers 
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have treated absorption [2l|, |22| and nonhnear dielectrics 
[H,!!!]. Knoester and Mukamel 0] and Huttner, Baum- 
berg, and Barnett [l^ start with the fundamental micro- 
scopic Hamiltonian and transform from coordinate space 
to wave-number space in the continuum approximation 
to derive the macroscopic Hopfield Hamiltonian [l^ . Re- 
cent work [13, m, [2^, [S^, [3l| includes macroscopic quan- 
tization of fields in magnetodielectric media, including 
left-handed negative-index materials. 

The spontaneous emission rate of an atom in a di- 
electric is readily derived by the macroscopic quantum 
electrodynamic theory. In 1976, Nienhuis and Alkemade 
used a macroscopic version of Fermi's golden rule to 
derive the dielectrically enhanced spontaneous emission 
rate with the macroscopic Ginzburg fields. Huttner, Bar- 
nett, and Loudon 32,], amongothers, have combined the 
macroscopic Fano-Hopfield [1, [1^ theory with Fermi's 
golden rule to derive the dielectric renormalization of the 
spontaneous emission rate of an impurity atom, while 
other studies begin with macroscopic Green's functions 
[131 or auxiliary fields [s^. For the most part, the col- 
lective attention is focused on the dielectric renormal- 
ization of the spontaneous emission rate. Knoester and 
Mukamel [3| obtained the dielectric effect on both the 
Lorentz redshift and the spontaneous decay rate by de- 
riving operator equations of motion from the Hopfield 
model. In this section, we derive equations of motion us- 
ing Weisskopf-Wigner theory applied to the macroscopic 
Hamiltonian in terms of the Ginzburg field operators '36| . 
A terni-by-term comparison with the generalized Bloch 
equations that were derived in the preceding section un- 
der the Lorentz local-field condition exposes an extraor- 
dinary degree of disagreement between two known and 
accepted treatments of the effect of a dielectric host on 
the electrodynamics of two-level atoms. 

The principal product of the macroscopic quantization 
theory is the medium-assisted field operator 
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where a| and a; are the macroscopic creation and destruc- 
tion operators for the field modes and mi is the frequency 
of the field in the mode I. Also, V is the quantization 
volume, GkjA is a unit vector in the direction of the po- 
larization, and A denotes the state of polarization. The 
spontaneous emission rate of an impurity atom in a di- 
electric can then be obtained by applying Fermi's golden 
rule. 
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to the effective interaction Hamiltonian Hmt = —fJ-a ■ E, 
where i labels the initial state and / denotes all avail- 
able final states. As is typically calculated, the dielectric 
renormalization of the vacuum spontaneous emission rate 
of an atom is found to be n [4, |5|, |6| due to the dielec- 
tric renormalization of the electric field operator (j3.ip by 



1/n, which is squared, and the D = density-of-states 
factor. 

Local-field effects of a dielectric are suppressed in the 
macroscopic quantization procedure and such local-field 
effects must be introduced phenomenologically. The 
paradigm that emerged from the propagation studies of 
Hopf and Scully [37t| and Bloembergen's work [lo| in non- 
linear optics is that the effect of a dielectric host is to 
multiply each occurrence of the dipole moment of a two- 
level atom by a local-field enhancement factor [1, [3, [13 ■ 
Using the Lorentz virtual cavity model of the local field, 
the spontaneous emission rate 
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for atoms in a dielectric scales as for large n. For an 
atom in a real cavity, the local-field enhancement factor 
is based on the Onsager model, and the modified spon- 
taneous emission rate 
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scales as n. A study of local-field effects by de Vries 
and Lagendijk [33] found that the Lorentz virtual-cavity 
model is appropriate if the atom goes into a crystal 
substitutionally but that the Onsager real-cavity model 
should be used for interstitial impurities. 

The macroscopic quantum electrodynamic theory can 
be used to derive additional consequences of the dielectric 
host for two-level atoms. Taking the field in a coherent 
state, the effective Hamiltonian is 
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(3.5) 

For species a, 0-3 is the inversion operator and crj_ are 
the raising and lowering operators for the j*'* atom, gj — 
{2nuji/hVy^'^^ai^j ■ ^kix) is the couphng between the 
atom at position rj and the radiation field, xj is a unit 
vector in the direction of the dipole moment at r^, u>a is 
the transition frequency, and /Iq is the matrix element of 
the transition dipole moment. 

Except for coefhcients of 1/n, and macroscopic field- 
mode operators, the effective Hamiltonian is the same as 
the microscopic Hamiltonian for identical two-level atoms 
in the vacuum. Equations of motion can then be de- 
rived in the same manner, which is the primary reason for 
adopting the macroscopic formalism. The formal integral 
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of the Heisenberg equation of motion for the field-mode 
operators is used to eliminate these operators from the 
remaining Heisenberg equations of motion. One obtains 
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factor of £ with each occurrence of the dipole moment_ 

based on the propagation studies of Hopf and Scully ^ ^ 

and Bloembergen's work in nonlinear optics. In this 
case, Eqs. p.8l) . become 
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and a similar equation of motion for the inversion op- 
erator. The procedure to evaluate these terms for two- 
level atoms in the vacuum is generally known and will 
be considered in detail in the next section. For now it is 
sufficient to note that the only differences from the vac- 
uum case are the coefficients of powers of n and the 
renormalization of the density-of-states in a dielectric in 
which the sum over modes is evaluated as 



dW 

dt 
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(3.9b) 

The ad hoc local-field correction gives the Bloembergen 
enhancement of the field that was derived classically in 
Eq. (|2.6p . although the redshift is over-corrected. More 
significantly, the Bowden-Lorentz redshift retains the ex- 
traneous factor of the refractive index due to the macro- 
scopically quantized fields. Then the equations of motion 
for two-level atoms in a dielectric host, Eqs. (I3.8|l and 
(13. 9p , that were derived using the macroscopically quan- 
tized fields are inconsistent with the generalized Bloch 
equations (|2.6p of the preceding section, and the appli- 
cation of quantized macroscopic fields to the electrody- 
namics of two-level atoms is contraindicated by classical 
Lorentz local-field theory. 
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in the mode continuum limit. Transforming to a rotating 
frame of reference and performing a local spatial average, 
one obtains the Bloch-likc equations of motion 
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where W = (cr3)sp and i?2i = {-i(J-)sp- 

Equations (j3.8p are the equations of motion for a two- 
level atom in a dielectric host medium that are derived 
using macroscopic quantum electrodynamics. The effects 
of the host appear as an enhancement of the decay rates 
and the Bowden-Lorentz redshift by a factor of the re- 
fractive index n, when compared to the vacuum ji — 1 
case. However, the equations of motion derived using 
the classical Lorentz local-field condition, Eqs. (|2.6p . dis- 
play an enhancement of both the Bowden-Lorentz red- 
shift and the field by a factor of £. 

The fact that local-field effects are suppressed in the 
macroscopic theory is well known. The accepted practice 
01 S [i3| is to phenomenologically associate a local-field 



IV. MICROSCOPIC QUANTUM 
ELECTRODYNAMICS 



In the preceding two sections, we derived generalized 
Bloch equations of motion for two-level atoms in a di- 
electric host by two well-known methods and obtained 
contradictory results. The most fundamental theoretical 
approach is to represent both the atoms and the dielec- 
tric microscopically 0, M, S M, 

and derive equations of motion from first principles. We 
show that the results of microscopic quantum electro- 
dynamics affirm the Lorentz local-field theory with re- 
spect to the field renormalization and the Lorentz red- 
shift. The macroscopic quantum electrodynamic theory 
produces different results for these effects and is therefore 
not valid. 

A dielectric host containing one or more two-level 
atoms is modeled quantum electromagnetically as a mix- 
ture of two species of atoms, a and b, embedded in the 
vacuum. To emphasize the symmetry of the local-field 
interaction, both species of atoms are initially treated as 
two-level systems. Species b is later taken in the har- 
monic oscillator limit that is associated with a relatively 
large detuning from resonance. Then, the total Hamilto- 
nian is comprised of the Hamiltonians for the free atoms, 
the free-space quantized radiation field, and the inter- 
action of the two-level systems with the free-space quan- 
tized electromagnetic field. The multipolar and minimal- 
coupling Hamiltonians are related by a canonical trans- 
formation and either can be used. However, due to 
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the canonical transformation, the circumstances of the 
rotating-wave approximation (RWA) are different for the 
two Hamiltonians [l^ . In the typical derivation of the di- 
electric susceptibility from the minimal-coupling Hamil- 
tonian, the RWA is invoked implicitly by replacing po- 
lariton eigenenergies with photon energies 0, M, Hsj . We 
take the direct route and use the multipolar Hamilto- 
nian in the RWA. Using a plane-wave expansion of the 
electromagnetic field 



E = V ^ («'^*""^ - ^■^■) ^k,A, (4.1) 



the multipolar RWA Hamiltonian is [11|, [4 
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dt 
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Bloch-like operator equations of motion are obtained by 
substituting the formal integral of the field-mode opera- 
tor equation (|4.3ap . 
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where a\ and ai are the creation and destruction opera- 
tors for the field modes and loi is the frequency of the field 
in the mode I with fc; = w;/c. For species a, is the in- 
version operator and (Tj_ are the raising and lowering op- 
erators for the J*'' atom, gl = (27rw//Sy)^/^/Lta(xj ■ ek,A) 
is the coupling between the atom at position rj and the 
radiation field, Xj is a unit vector in the direction of 
the dipole moment at rj, uja is the transition frequency, 
and fia is the dipole moment. For species b, <;^, <;jj, /i", 
r„, Wfc, and fi^^ perform the same functions. Also, V is 
the quantization volume, ek,A is the polarization vector, 
and A denotes the state of polarization. The polariza- 
tion indices on the variables have been suppressed for 
clarity. In the two-level approximation, the transitions 
s — Am e (—1,0,+!) are treated separately and the op- 
erators need not carry a specific value for the magnetic 
sublevel [i^ . 

Equations of motion for the material and field-mode 
operators are developed in a straightforward manner 
from the Hamiltonian using the Heisenberg equation. We 
have 
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into the material operator equations of motion (|4.3bp . 
(|4.3cl) . (I4.3dp . and (|4.3ep [i^]. We transform operator 
variables to different rotating frames of reference in which 
a-'_ = cr^e*'^"* and = c^^e*"''* are slowly varying quan- 
tities. Performing the indicated substitution into Eq. 
(|4.3bp produces 



da^ 

^=E3f-i(0e''''-'a.(0)e 

IX 



(A 



"'0 



IX 



^h\'*q1{t')e'^' 



(r,-r.) 



(4.5) 



in normal ordering with Vji — rj — Vi. Likewise, one 
obtains 
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Then 



from Eq. (j4.3dp . Equations of motion for the inversion 
operators are obtained from Eqs. ()4.3cp and (|4.3e|) in a 
similar fashion. 

The field that drives an atom, Eq. (|4.4p . consists of 
the vacuum field, the self-field, and the reaction field and 
we can identify the terms on the right-hand side of Eqs. 
(|4.5|) and (|4.6p with fluctuations due to the vacuum field, 
spontaneous decay from the self-field, near-dipole-dipole 
interactions between same-species atoms associated with 
the reaction field, and near-dipole-dipole interactions of 
an atom with the atoms of the other species, also asso- 
ciated with the reaction field. The usual procedure is to 
limit consideration to only the spontaneous decay rate 
of a single impurity atom by dropping the fluctuations 
and the single-species interactions for both the bath and 
impurity atoms. These terms are retained here because 
they contain significant information about local-field ef- 
fects in dielectrics. For concreteness, we take species b 
to be the bath atoms and species a to be the two-level 
impurity atoms. 
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t plying the typical Weisskopf-Wigner procedure |49l . 
|51| in the mode continuum limit, one obtains 
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where 
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is the spontaneous decay rate. For an atom of species 
6, initially in the excited state, 7h is the spontaneous 
emission rate Fq into the vacuum. 

The pairwise interaction of atoms is carried in the re- 
maining m ^ n part of the summation. In the Milonni- 
Knight [46l| model of the interaction of two identical two- 
level atoms, the strength of the interaction depends on 
the separation distance and the magnetic sublevel tran- 
sition. Then 



V. NEAR-DIPOLE-DIPOLE INTERACTION 

The near-dipole-dipole interaction is the basic mech- 
anism of the action of the local field. The Weisskopf- 
Wigner-based procedure to evaluate the dipole-dipole 
interaction for a dense collection of identical two-level 
atoms was developed by Ben-Aryeh, Bowden and En- 
glund [Till, with corrections by Benedict, Malyshev, Tri- 
fonov, and Zaitsev 48], to investigate single-species in- 
trinsic optical bistability. The results apply to both 
species of two-level atoms, individually, but we work with 
species b in order to maintain consistent notation when 
we take the harmonic oscillator limit of a two-level atom 
and derive the Lorentz local-field correction in a dielec- 
tric. 

We consider a dense collection of identical two-level 
atoms of species b in which the atoms are evenly dis- 
tributed in the vacuum with a number density TVf,. The 
same-species interaction 



Fi(i?) = e*^ | -- + — + — 
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for Am = ±1 transitions and 
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for Am = transitions, where 
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R — ki,r,i„i — LJi,rnm/c, kb — |kb|, Vnm — !"„ , and 
fnm — \^nm\- Performing the summation over the mag- 
netic sublevels, the pairwise dipole-dipole interaction can 
be written as 

It" = ^3 E Bnm^"^{t - r„,„/c), (5.8) 



where 



X ^/i J" )e'''''^'""~'"'" ^ (5-1) 

vas 

can be extracted from Eq. (|4.6p . 

The self-interaction of the n*'' atom with its own re- 
action field is characterized by the term r™ — r„ in the 
interaction as a consequence of the relation <j|'(t)<?" (i) = 
— <J"(t) between Pauli spin operators for the same atom. 



-I- (X,„ • n„„i)(x„ • n.nm)F2{hrnm) (5.9) 

and hnm = Ynm/rnm is a unit vector in the direction of 
Trim = !■« - Tm [13, El ■ Further, Bnm incorporates a 
view factor to account for the arrangement of the dipoles 
in the volume. 

The atoms are evenly distributed with a number den- 
sity Ni,. For the n*'* atom, the single-species dipole- 
dipole interaction is obtained in a summation over all 
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other atoms of species b. In the region near r„ , the inter- 
action is evaluated by taking the location of dipoles as 
discrete, while the continuum approximation is applied 
elsewhere. Then 



V-Vs 



where S is the radius of a small spherical volume Vg , larger 
than a cubic wavelength, about the point r„. For cu- 
bic symmetry, the field generated by the localized atoms 
Tj, in the virtual cavity is zero at the center 

ullilEi- 

The atom n is located at the origin of a cylindrical 
volume of thickness L and radius Rq. The near-dipole- 
dipole interaction is obtained by evaluating the integral 





fL/2 j- 


/ d4 


/ dz 


Jo 
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Fi + cos'^ (f>F2^, (5.11) 



excluding a volume {A/3)Trd^ about the origin from the 
range of integration, resulting in 

J,, ^ Z^Illdl^e^k,s^n^t){^^{t~\r\/c))s^ (5.12) 
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In the limit S 0, the ncar-dipole-dipole interaction 

It" = -iiyt^Sit)^- (5.13) 
remains finite. Here 

An 



(5.14) 



is the strength of the near-dipole-dipole interaction and 
■f- = {^-[t — |r|/c))sp represents a spatially averaged 
quantity. The details of this calculation can be found 
in the articles by Ben-Aryeh, Bowden, and Englund 14| 
and by Benedict, Malyshev, Trifonov and Zaitsev [48 1. 

The atoms of species b can be treated as harmonic 
oscillators if all excitation frequencies are far from res- 
onance with LOb- In this limit the atom essentially re- 
mains in the ground state such that ^3 — > —1. Then 
the near dipole-dipole interaction reduces to the Lorentz 
local-field correction, shifting the resonance frequency by 
47riVb|^fcp/(3S.). The microscopic result is in full agree- 
ment with the classical Lorentz local-field correction and 
has been experimentally validated [s^ . Issj by selective 
reflection of Rb from a sapphire window. 



Finally, all of the results of this section can be applied 
to the other species of atom. Repeating for species a 
yields 



If « -iVa(Ti{t){a^{t-\v\/c)),p = -iva4{i)^-- (5-16) 



Bq^{t-\r\/c)d^Y, (5.10) In addition. 
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la = 2/3, 



and 



At: 



(5.17) 



(5.18) 



are defined for later use. 



VI. INTERSPECIES INTERACTION 

The effect of the interspecies near-dipole-dipole in- 
teraction can also be evaluated microscopically using 
Weisskopf-Wigner theory. A single rotating frame of ref- 
erence is used for both species of atoms by making the 
transformation ^" — q^^Q-^i^b-^a)t ^ Applying the results 
of the preceding section, we have 



~dt 



(6.1) 

in the harmonic oscillator limit <^3 (i) —1, where the 
fluctuating field 



= ^(2^cc>,ft/y)i/2(x„ . efe)e'''--"az(0)e 



-i(aj;-!^b)t 



l\ 



is associated with the spontaneous decay rate by the 
Kramers-Kronig relations. 

The interspecies interaction in Eq. (|6.ip describes how 
a specific host atom n interacts pairwise with each of the 
impurity atoms. The formal integral of the equation of 
motion of the host atoms, Eq. (|6.ip . is 

Jo ri Jq 

xJ2K f di"e-^('"'-'"'')(*'-*")^5raL(i")e^*''-'"% 
;a -^0 is 

(6.2) 
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where a = LUf, — uja ~ h'b — ^76/2. Substituting Eq. (|6.2 
into Eq. (|4.5p . one obtains 



atoms in the continuum hmit. For nonidentical atoms, 
the pairwise interaction goes as [13, 



da' 



where 



^"0-3/ - ii^aO-gCT- - ^o-i + /2, (6.3) 



xe 



t" 



dt" 



+ 7;VKfhNa Ba-{t" -\r\/c)d^r. (6.6) 



Then, Eq. (|6.6p is evaluated as in Sec. V to obtain [ij 



rdd^ 



3h 



NafilM^-it" -\r\/c)),p 



(6.4) 



The term containing the fluctuations has been dropped 
from consideration because the procedures presented here 
apply only to slowly varying quantities and because there 
will be no contribution from the random fluctuations af- 
ter averaging. 

The direct dipole-dipole interactions between impu- 
rity atoms was derived in Section V. Equation (|6.4p con- 
tains two such dipole-dipole interactions between non- 
identical atoms that are integrated over the different 
subspaces corresponding to (i) bath atoms and (ii) two- 
level impurity atoms. Due to the relation a'^{t)a-'_{t) — 

—a-L{t) between Pauli spin operators, the term i = j is 
the special case that is associated with the renormaliza- 
tion of the spontaneous decay rate. This separates the 
interaction I2 = 12'^^ + into for the case i = j 
and 12'^'^ for the summation over the rest of the impurity 
atoms. The two parts of the interaction will be consid- 
ered separately. 



A. Dielectric Mediated Dipole— Dipole Interaction 

The n*^ atom of the dielectric interacts pairwise with 
every impurity atom. The term 



rdd-' 



I'X' -^0 



^^///g-i(w,,-l^a)(t"-t"') 



it" 



„ik,/-(r„-r. 



(6.5) 



extracted from Eq. (|6.4p . can be evaluated in the same 
manner as in Section V, except that the atoms are of dif- 
ferent species. The summation represents the effect of all 
the impurity atoms on a single atom of the host material. 
The sum over the i ^ j impurity atoms is performed (i) 
in the near region by the discrete summation over the im- 
purity atoms and (ii) elsewhere by treating the impurity 



47ri ^ _ 



(6.7) 



The quantity ct_ = {^-{t" ^ kl/c))sp is slowly varying in 
time and the temporal integral 

ir'it') = dt"e-(*'-*") 7V„M>.^- (6.8) 

can be performed in the adiabatic-foUowing approxima- 
tion. Repeatedly integrating Eq. (|6.8p by parts ^56*1, the 
series can be truncated at the first term in the expansion 
if the time rate of change of (t_ is much smaller than a(7_ 
yielding 



ir it') 



1 \ /-47ri 



3^ 



Na^J.l^J.bcr- 



(6.9) 



The remaining part of the I2 integration is another 
interspecies dipole-dipole interaction. In this case, the 
summation imparts the effect of all the atoms of the host 
dielectric, modified by interspecies interaction with the 
impurity atoms, on the j'*'' impurity atom. Combining 
terms, 

l2''{t') = ^mfilfia-^Nafiytaia. (6.10) 
6h a An 



becomes 



l2\t') = -^Xbi^aaia., (6.11) 



where Xb is the linear susceptibility of species b. Adding 
the direct near-dipole-dipole interaction from Eq. ()6.3p . 
we obtain 



-f 2 

-i — - — i/qCTs {t)a^ (t) . 



(6.12) 



Comparison of Eq. (|6.12p with the single species dipole- 
dipole interaction, Eq. (|5.16p . shows that the effect of the 
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dielectric host is to enhance the interaction by a factor 
of £ — {-n? + 2)/3. Taking the local spatial average, W = 



(0-3 )sp and R21 
redshift 



(— i(T_)sp, one finds that the Lorentz 



+ 2 477 

3 m 



(6.13) 



is consistent with the Lorentz local-field calculation, Eq. 
(|2.6ap . while the macroscopic quantum electrodynamic 
result, Eq. (|3.9ap . is not. 



B. Dielectric-Enhanced Spontaneous Decay Rate 

Most of the elements of the microscopic theory of the 
spontaneous decay rate of an atom in a dielectric are 
common to the treatment of the dipole-dipole interac- 
tion. In order to show this clearly, we consider an equiv- 
alent derivation of the dielectric mediated dipole-dipole 
interaction. Performing the temporal integrations first, 
the interspecies interaction (|6.4p can be written as 



six 



k,,-(r„-ri 



'V\' 



(6.14) 

Applying the Milonni-Knight interaction with a view fac- 
tor results in 

h = - — ^P'^Pbcri ^\[{S^j ■ n„j)(x„ • n„j)]F2 (fcar„j) 



[(xj • x„) - (xj • n„j)(x„ • n„j)] Fi{karnj)^ 



integration over the subspace of two-level atoms. Like- 
wise, the summation over the magnetic sublevels invokes 
5ss'- Then, 



+ [(% ■ X„) - (Xj ■ n„j)(x„ • nnj)f i^f (fcafnj)|- (6.16) 

The microscopic treatments of the spontaneous decay 
rate in a dielectric [lllj 13] are missing elements of Eq. 
(j6.16p and can neither affirm nor contradict the macro- 
scopic theory of quantum electrodynamics. Because the 
dielectric renormalization of the spontaneous decay rate 
does not have a classical local-field condition-based ana- 
log, we do not consider it further. 



VII. DIELECTRIC-ENHANCED FIELD 

The dielectric has an effect on an applied electromag- 
netic field that can also be evaluated microscopically. 
Taking the field in a coherent state, the partial Hamilto- 
nian is 



ih X - 



n 

(7.1) 

where cUp is the nominal frequency of the field and fla = 
pLaS/h and rj;, = fif,£/h are Rabi frequencies. The total 
Hamiltonian is now comprised of the Hamiltonians (j7.ip 
and (|4.2p . Developing Heisenberg equations of motion 
and eliminating the field-mode operators and the dielec- 
tric operators results in the appearance of 



^ |[(xi • n„i)(x„ • n.ni)]F2{kar„i) 



nslX 



+ [(x,; • X„) - (x,; • n„i)(x„ • Uni)] Fi (fcar^) jcrL . (6.15) 

Converting the sums, excluding i — j, to integrals and 
integrating over the subspace of two-level atoms and then 
over the subspace of oscillators is equivalent to the deriva- 
tion of the dielectric mediated dipole-dipole interaction 
that was presented in the preceding subsection. 

The renormalization of the spontaneous decay rate of 
a dielectric-embedded two-level atom is derived from the 
interspecies interaction (j6.4p in the same fashion by tak- 
ing the target atom to be the same as the source atom. 
The summation over the two-level atoms is evaluated 
with the use of the delta- function (J,,, rather than the 



as an addition to Eq. (|6.3p . Equation (|7.2p contains the 
same type of interaction that was evaluated in the previ- 
ous section. Performing the adiabatic-foUowing approxi- 
mation and the sum over polarizations, bath atoms, and 
magnetic sublevels in the mode continuum limit, we ob- 
tain 



n'^ + 2fiar- j 

cai. 

3 2h 3 



(7.3) 



The electromagnetic field is enhanced by the same factor 
of £ as the reaction field. 
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VIII. OPTICAL BLOCK EQUATIONS FOR 
EMBEDDED ATOMS 



The macroscopic optical Bloch equations can be de- 
rived from the quantum electrodynamic equations of mo- 
tion in the Hmit of large numbers. Combining Eqs. (|6.3p . 
(|6.1ip . and (|7.2I1 . one obtains 



-liiyaCr'a^+i^Se-'^^-'^-'^'a^:,, (8.1) 
dt ^ 2h 3' V ; 



neglecting the Gaussian noise source with zero mean and 
absorption. The equation of motion for the inversion 
operator 



dal 
dt 



H.c. 



is derived in a similar manner. Optical Bloch equations of 
motion are obtained by transforming to a frame rotating 
at the frequency of the field and taking a local-spatial 
average, as in Sec. III. We compare the optical Bloch 
equations 



dR2i 
dt 



47r 
3^ 



2h 



dW 
'df 



[fiats* R2I - fiae£R*2i] 



.3a) 



..3b) 



that were derived from first principles, to the Lorentz 
local-field-based equations (|2.6p . Based on a favorable 
comparison of the local-field enhancement of the Lorentz 
redshift and the Rabi frequency with the classically de- 
rived result, we can reasonably assert that the micro- 
scopic theory, unlike the macroscopic quantum electro- 
dynamic theory, satisfies the correspondence principle. 

The microscopic theory allows us to consider the more 
general case of of a complex local-field enhancement fac- 
tor. Separating the real and imaginary parts of £, the 
optical Bloch equations can be written as 



dR2i 
dt 



2h 



OJp-UJa- ^NfllirW ] R21 



££W ~^Nnle^WR2i, 



(8.4a) 



dW 
'df 



[flat£*R2l- ^la^£R*2l]-^^^ya\R21? (8.4b) 



with i — + iti. The microscopic theory justifies the 
use of a complex refractive index in the classical Lorentz 
local- field condition. Then Eqs. (|8.4p . with phenomeno- 
logical damping, can be derived by substituting the po- 
larization (|2.5p with complex n into the generalized Bloch 
equations (|2.4p . The imaginary part of the Lorentz red- 
shift, derived in this manner, was found to be associated 



with an intrinsic cooperative decay for two-level atoms 
in an absorptive host @ . This result is confirmed by the 
microscopic theory. 

The optical Bloch equations (|8.4p for dielectric- 
embedded two-level atoms are derived from the mi- 
croscopic description of quantum electrodynamics using 
vacuum-based fields that are known to satisfy the equal- 
time commutation relations. Because the field-mode op- 
erators have been eliminated, the equal-time commuta- 
tion relations cannot be discussed in the context of the 
optical Bloch equations (|8.4p or Heisenberg equations 
(|8.ip and (|8.2p . Instead, the optical Bloch equations, gen- 
eralized for a dielectric host, must demonstrate conserva- 
tion of probabiHty. The total population is W^-|-4|i?2iP- 
Direct substitution from Eqs. (|8.4p shows that the tem- 
poral derivative of this quantity is nil, as required, in the 
limit that absorption by the atoms and the host dielectric 
can be neglected. 



IX. SUMMARY 

The interest in the dielectric renormalization of the 
spontaneous emission rate of an atom embedded in a di- 
electric material has obscured the inconsistencies in the 
macroscopic theory of quantum electrodynamics. The 
dielectric renormalization of the Lorentz redshift and the 
Rabi frequency, but not the spontaneous decay rate, can 
be derived using the classical Lorentz local-field con- 
dition providing an independent check on the validity 
of macroscopic quantum electrodynamics. The optical 
Bloch equations for a dense collection of two-level atoms 
in a dielectric host medium were derived using the clas- 
sical Lorentz local-field condition 

47r_ 



E, 



E 



3 



in Sec. II, while in Sec. Ill, a different set of optical Bloch 
equations were derived using the macroscopic quantum 
electrodynamic theory. Both derivations are short and 
uncomplicated and are based on well-established phys- 
ical theories, yet lead to contradictory results for the 
Lorentz redshift and the Rabi frequency. If we assume 
the validity of the Lorentz local-field condition, then the 
macroscopic procedure is proven to be incorrect. Con- 
versely, the validity of the macroscopic quantum elec- 
trodynamic theory would imply that the Lorentz local- 
field condition is incorrect. One deciding factor is that 
the Lorentz local-field correction has been validated ex- 
perimentally [s^, ^535 , while the experimental record for 
the macroscopic quantum theory has been inconclusive. 
We applied the more fundamental microscopic theory 
of quantum electrodynamics to the same problem and 
demonstrated complete agreement with classical theory. 
The differences in the Rabi frequencies can be reconciled 
with a phenomenological local-field factor applied in the 
macroscopic case, providing the virtual-cavity model is 
used. However, no such facile reconciliation can be pro- 
vided for the Lorentz redshift. We conclude that both 
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the correspondence principle and microscopic quantum tum electrodynamic theory, 
electrodynamics are violated by the macroscopic quan- 
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